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We analyze the finite size corrections to entanglement in quantum critical systems. By using 
conformal symmetry and density functional theory, we discuss the structure of the finite size con- 
tributions to a general measure of ground state entanglement, which are ruled by the central charge 
of the underlying conformal field theory. More generally, we show that all conformal towers formed 
by an infinite number of excited states (as the size of the system L — > oo) exhibit a unique pattern 
of entanglement, which differ only at leading order {1/ Vf . In this case, entanglement is also shown 
to obey a universal structure, given by the anomalous dimensions of the primary operators of the 
theory. As an illustration, we discuss the behavior of pairwise entanglement for the eigenspectrum 
of the spin- 1/2 XXZ chain with an arbitrary length L for both periodic and twisted boundary 
conditions. 
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I. INTRODUCTION 

In recent years, the observation that entanglement 
may play an irnportant role at a quantum phase tran- 
sition [l|, 0, [1, |j] has motivated intensive research on 
the characterization of critical phenomena via quantum 
information concepts. In this direction, conformal invari- 
ance has brought valuable information about the behav- 
ior of block entanglement, as measured by the von Neu- 
mann entropy, in critical models. Indeed, conformal field 
theory (CFT) has been used as a powerful tool to deter- 
mine universal properties of entanglement. Remarkably, 
it was shown that the entanglement entropy obeys a uni- 
versal logarithmic scaling law for one-dimensional criti- 
cal models both at zero and finite temperatures [1, d, 0| j 
which is governed by the central charge of the associ- 
ated CFT. Moreover, corrections to the entanglement 
entropy due to finite size effects have also been consid- 
ered [a, Q for periodic and open boundary conditions. 
Together with approximative methods such as renormal- 
ization group (see, e.g. Refs. lol. [lOl . [ill . [T^. [l3t) and den- 
sity functional theory (DFT) [ij, CFT has been settled 
as one of the most promising approaches for investigat- 
ing the behavior of entanglement in many-body quantum 
critical systems. 

In this work, we will exploit in a new perspective the 
impact of CFT methods for the evaluation of entangle- 
ment at criticality. More specifically, our approach will 
be based on the statement that finite size corrections to 
the ground state expectation values of arbitrary observ- 
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ables are ruled by conformal invariance. This conclusion 
is indeed a consequence of two results: (1) Finite size cor- 
rections to the energy spectrum of a critical theory are 
determined by conformal invariance [H,[lB,[l3|; (2) DFT 
techniques imply that, under certain conditions discussed 
below, general observables can be evaluated as a function 
of the first derivative of the ground state energy with re- 
spect to a Hamiltonian coupling parameter flS, . We 
then simultaneously apply these two results to obtain 
the finite size corrections to ground state entanglement 
in critical models. As a by-product, conformal invari- 
ance determines the structure of entanglement in the 
presence of extra symmetries for certain higher energy 
states, which are the lowest energy states in each sym- 
metrically decoupled subspace of the Hilbert space. For 
instance, if the Hamiltonian is translationally invariant 
and has a U{1) symmetry due to its commutation with 
the magnetization operator, we can split out the Hilbert 
space in sectors of fixed momentum and magnetization. 
More generally, we will also show that all conformal tow- 
ers formed by an infinite number of excited states (as the 
size of the system L — > oo) exhibit a unique pattern of 
entanglement, which differ only at leading order (1/i)^. 
This will be based on a generalization of the HK theo- 
rem for individual states belonging to conformal towers 
of critical systems. Finite size corrections to entangle- 
ment in these excited states will obey a universal struc- 
ture, given by the anomalous dimensions of the primary 
operators of the theory. 

Since our approach is applicable for any entanglement 
measure, it allows in particular for the investigation of 
the universality properties^of pairwise entanglement mea- 
sures, e.g., concurrence [20] and negativity [21]. For pair- 
wise measures, criticality was first noticed through a di- 
vergence in the derivative of entanglement, signaling a 



2 



second-order phase transition For first-order phase 
transitions, jumps in entanglement itself indicates quan- 
tum critical points [13, [l^- A general explanation for this 
distinct usual behavior of first-order and second-order 
phase transitions has been provided in Refs. [l^, (for 
an explicit discussion of examples which do not obey this 
expected behavior, see Ref. [25{). From the point of view 
of CFT, we will be able to explicitly work out the finite 
size corrections to pairwise entanglement measures and 
show how these corrections involve universal quantities, 
such as the central charge or the anomalous dimension of 
primary operators associated with the CFT. As an illus- 
tration, we will consider the spin-1/2 XXZ chain, where 
an analytical expression, valid up to o(L^^), will be pro- 
vided for the negativity of nearest neighboring spins as a 
function of the size L of the chain. 



II. ENERGY SPECTRUM AND FINITE SIZE 
EFFECTS IN CRITICAL QUANTUM SYSTEMS 

Let us consider a critical theory in a strip of finite 
width L with periodic boundary conditions. The transfer 
matrix of the theory is written as T = exp{aH), where 
a denotes the lattice spacing and H is the Hamiltonian. 
Then, for large L, the ground state energy density e(L) = 
En{L)/L of H is provided by conformal invariance jlSl . 
[l6j . reading 

e{L)=e^-^L-^+o{L~^), (1) 
6 

where Soo is the energy density in the limit L — > cx) and 
o{L~'^) denotes terms of any order higher than . In 
Eq. ([T]), c is the central charge of the Virasoro algebra 
(the conformal anomaly) and the parameter ^ must be 
fixed in such a way that the equations of motion of the 
theory are conformally invariant [26j . The structure of 
the higher energy states is determined by the primary 
operators of the theory For each operator Oa with 
anomalous dimension Xa, there corresponds a tower of 
states with energy densities e"j, {L) given by 

el^,{L) = e{L) + 2ti i{x^ + j + f)L-^ + o{L-^), (2) 

where j,j' = 0,1,... are indices labelling the tower of 
states associated with the anomalous dimensions Xa- 
Higher-order corrections to Eqs. ^ and ^ as well as 
convenient generalizations for more general boundary 
conditions, e .g.. twisted boundary conditions, may also 
be obtained 

III. HOHENBERG-KOHN THEOREM AND 
EXPECTATION VALUES OF OBSERVABLES 

Let us turn now to the discussion on how DFT can 
be allied with conformal invariance to extract informa- 
tion about expectation values of observables from the 



energy spectrum. DFT [2^ [s^l is originally based on the 
Hohenberg-Kohn (HK) theorem [2^ which, for a many- 
electron system, establishes that the dependence of the 
physical quantities on the external potential v{r) can be 
replaced by a dependence on the particle density n(r). 
The HK theorem can be extended for the context of a 
generic quantum Hamiltonian H on a lattice (see, e.g., 
Refs. [l8, 19]). In order to be specific, let us consider a 
quantum spin chain of size L governed by the Hamilto- 
nian 

L 

H = Ha + \Y,A,, (3) 

1=1 

where A is a control parameter associated with the Hermi- 
tian operators Ai which act on the site i, e.g., an observ- 
able relevant to driving a quantum phase transition. Let 
us take, for simplicity, a translationally invariant chain 
(e.g., by assuming periodic boundary conditions). Then, 
by taking the expectation value of Eq. ([3]), we obtain 

{H)^{Ha)+\L{A), (4) 

where {A) = (yl^) = (Aj) {Vi-J) due to translation sym- 
metry. Therefore 

e = eo + A(yl>, (5) 

where e = {H) /L and Eq {Ho) / L are the energy densi- 
ties associated with H and Hq, respectively. For a general 
Hamiltonian such as given in Eq. ([31), the HK theorem 
can be generalized to the statement that there is a du- 
ality (in the sense of a Legendre transform) between the 
expectation value {A) (corresponding to n(r)) and the 
control parameter A (corresponding to v{y)) [i1,[i^. In 
order to specify the conditions supporting this duality 
let us separately consider the cases of nondegenerate and 
degenerate Hamiltonians. 

A. Nondegenerate case 

Let A and A' be two fixed values of the coupling pa- 
rameter in Eq. ([3]), which correspond to nondegenerate 
ground states given by 1-0) and respectively. We as- 
sume that, for A ^ A', we have that \ip) ^ with 
a a complex phase. This assumption means that differ- 
ent values of the coupling parameter are associated with 
distinct ground states. It reflects the requirement of the 
uniqueness of the potential (see, e.g., Ref. [3l|). A gen- 
eral condition to ensure the uniqueness of the potential 
for Hamiltonian ([3]) will be derived below. Then, by as- 
suming a unique potential and taking two different cou- 
plings A and A', the Rayleigh-Ritz variational principle 
allows us to write 

{mm < {i^'iHm = {i^'W) + (a - ao L{Ay , (6) 

where {A)' = {tp'\A\il!') and H and H' are the Hamilto- 
nians associted with A and A', respectively. Therefore 

e-e' < (A-A')(^)'. (7) 
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Analogously, application of the variational principle for 
the ground state results into 



e' -e< (A' -A) {A). 
By adding Eqs. ^ and fS]) wc obtain 

{A)'^{A). 



(8) 



(9) 



Eq. ([9]) expresses the HK theorem for nondegenerate 
ground states, stating that distinct densities are asso- 
ciated with distinct potentials. In other words, we can 
establish the map 

A ^ IV-) ^ {A) = . (10) 



B. Degenerate case 

In order to establish the HK theorem for degenerate 
ground states, let us consider two fixed values of the cou- 
pling constant, each of them associated with arbitrarily 
degenerate ground states: 

A < — > q — degenerate ground states: • ■ • j IV'g)} ■, 

A' < — > g' — degenerate ground states: l^/iq/)} . 

Considering that any of the ground states are equally 
likely, we can describe them by the uniformly distributed 
density matrices 



(11) 



The requirement of uniqueness of the potential yields in 
the degenerate case the condition that p and p' are dis- 
tinct. Applying the variational principle, we obtain 

Tr {pH) < Tr {p H) = Tr {p'H') + {X - A') , (12) 

where, here, (A)' = Tr{p'A). Eq. ^ implies that 
£ — e' < (A — A') {A)'. Therefore, as before, we use the 
complementary equation e' — e < (A' — A) (A) and obtain 
(A)' 7^ (A). The HK map in this case can be written as 



Aj, as given in Eq. do not have common eigen- 
states. Sufficiency: Suppose that two distinct couplings 
A and A' yield the same eigenstate of H 



Then, from Eqs. and we obtain 



(14) 



(15) 



(16) 



Therefore, in this case, 1-0) is also an eigenstate of J^i 
(as well as an eigenstate of Hq). Hence, the condition 
that Hq and J2i ^'^^ have common eigenstates is 

sufficient for ensuring the uniqueness of the potential. 
Necessity: Let us suppose that Hq and J2i have a 
common eigenstate 



(17) 
(18) 



Then we obtain that {Hq + A X), A,)\ip) = [Eq + Aa)|V'). 
Hence, by varying A, we only change the eigenvalue (keep- 
ing the same eigenstate) , which means that distinct cou- 
plings will lead to the same eigenstate of H . Therefore, 
the condition that Hq and Ai do not exhibit a com- 
mon eigenstate is also necessary for the uniqueness of 
the potential. In conclusion, the sufficient and neces- 
sary condition for the uniqueness of the potential can be 
translated by the noncommutation relation 



(19) 



Naturally, we disregard in Eq. (|19p the rather nonusual 
situation where Hq and Ai are noncommuting observ- 
ables, but [-f^o, X)i ^i]IV') results in a vanishing quantum 
state. 



{A)=TT{pA) 



(13) 



D. HK theorem for conformal towers in quantum 
critical models 



C. Uniqueness of the potential 

As discussed above, the condition for the uniqueness 
of the potential, which is fundamental for the derivation 
of the HK theorem, is defined by the requirement that 
different values of the coupling parameter A are associ- 
ated with distinct ground states of the Hamiltonian H . 
Here we will show that a necessary and sufficient con- 
dition for which different values of X are associated with 
distinct eigenstates of H is that the operators Hq and 



Since the HK theorem is based on a variational princi- 
ple, we cannot guarantee that the expectation values of 
the observables in individual excited states are in general 
a function of the derivative of the energy of the excited 
state. Naturally, as previously mentioned in Sec. HI the 
HK theorem can be applicable in the presence of symme- 
tries to excited states that are the minimum energy states 
in a given symmetric subspace of Hilbert space. In this 
work, we show that, under certain conditions, the HK 
theorem can also be extended for all the individual states 
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of conformal towers in quantum critical models. We be- 
gin by supposing a periodic chain governed by a Hamilto- 
nian given by Eq. ([3]) which is conformally invariant in a 
critical interval < A < • Moreover we will assume 



the condition pQ]) for the uniqueness of the potential. Let 
us denote by {|V'^j/.rf(A))} the set of cigcnstates associ- 
ated with the energy e"j, (A), with d— 1, . . . ,D labeUing 
the 15-fold degeneracy (see Eq. ([2|)). We take the system 
in a uniformly distributed density matrix 



(20) 



d=l 



Our aim is to show that the potential A uniquely specifies 
the density 

mr-T^[f^,rm^^. (21) 



Therefore, the derivative 9e"y/9A must be a monotonic 
function of A. In order for this to occur, it is sufficient 
that: (i) de'^j,/dX is continuous in the interval < 
A < Xc2 and (ii) 9^£"^,/9A^ ^ 0. Condition (i) is usually 
achieved for a smooth (well-behaved) energy. Concerning 
condition (ii), let us take the derivative of Eq. which 
yields 



2_Q 
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The first term in the r.h.s. of Eq. ([22]) concerns the sec- 
ond derivative of the ground state energy with respect 
to A. We can show that this term is strictly negative. 
Indeed, from Eqs. ([7]) and ([5]), which hold for both de- 
generate and non-degenerate ground states, we obtain 



(23) 



where (A) denotes the expectation value of A taken in the 
ground state. Concerning the second term in the r.h.s. 
of Eq. (P^ . it is negligible for large L. Consequently, we 
can write 



dX^ 



U < (large L). (24) 



Hence, d'^e" j, /dX^ is non- vanishing and then the deriva- 
tive de"j,/dX is monotonically related to A. Therefore, 
a D-fold degenerate (up to order L~^) eigenlevel given 
by a, j, and j' defines a density matrix p"j, that can be 
taken either as a function of A or {A)'jj, . This provides 
an extension of the HK theorem for arbitrary individual 
eigenstates belonging to conformal towers in quantum 
critical models. 



specially useful for the investigation of entanglement in 
the ground state of quantum systems undergoing quan- 
tum phase transitions [l^. In particular, the dependence 
of an arbitrary entanglement measure M on the param- 
eter A can be replaced by the dependence on the ground 
state expectation value {A) [l^, which means that 

M = M(A) = Mi{A)) = M(|i), (25) 

where the Hellmann-Feynman theorem (3^ . [ssj has been 
used in the last equality above. As discussed in the last 
Section, in the case of critical models, the HK theorem 
can also be applied to any state of conformal towers, 
which allows us to write the entanglement of such states 
as 



Eq, 



M-r-M^,{X} = Aq,{{A)l,). 
can be rewritten by observing that 



(26) 



(27) 

and performing a 



By inserting Eq. (P7)) into Eq. 
series expansion, we obtain 



dX 



(28) 



This means that the entanglement corresponding to all 
the eigenstates in the conformal towers are nearly the 
same as that of the ground state, with corrections of or- 
der . Moreover, it follows that finite size effects for 
an arbitrary measure of entanglement are ruled by con- 
formal invariance. 

In order to evaluate entanglement at a point A = Ac, 
we should be able to perform a derivative of the energy 
with respect to A taken at Ac- Therefore, note that we 
will have that Eqs. H]) and ^ are the starting point to 
determine finite size corrections to entanglement if and 
only if the theory is critical in an interval around Ac. 
For the case of a single critical point (e.g., Ising spin-1/2 
chain in a transverse field) instead of a critical region 
(e.g, XXZ spin-1/2 chain in the anisotropy interval — 1 < 
A < 1), more general expressions for the energy should 
be used, which take into account a mass spectrum (see, 
e.g., Ref. (H). 



V. ENTANGLEMENT IN THE FINITE SIZE 
SPIN-1/2 XXZ CHAIN 



IV. FINITE SIZE CORRECTIONS TO 
ENTANGLEMENT IN CONFORMAL 
INVARIANT MODELS 

The HK theorem implies a duality between the poten- 
tial A and the density {A). This behavior was revealed 



As an illustration of the previous results, let us con- 
sider the spin-1/2 XXZ chain, whose Hamiltonian is given 

by 

Hxxz = - 2 E («+i + ^'^'+1 + ' (29) 

i=l 
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where periodic boundary conditions (PBC) are assumed. 
We will set the energy scale such that J = 1. Entan- 
glement for spin pairs can be quantified by the negativ- 
ity [2l|, which is defined by 



Af{p'^) = 2 max(0, - min(/ijf )), 



(30) 



where /x^ are the eigenvalues of the partial trans- 
pose p^i''^^ of the density operator p*-', defined as 
{a(3\p^^ 17,5) = {'^P\p\a5). For the XXZ model, U{1) 
invariance ( [ff, ct*] = 0) and translation invariance 
ensure that the reduced density matrix for spins i and j 
reads 



( \ 

b'^ 


\ d' ^ ) 



(31) 



where 



-(1 + 2G, 

\ (1 - ge: 



= l[iGt + G'^^+^iG%-G^^], (32) 

where Gz — {a],) is the magnetization density (computed 
for any site i) and G^^ = (c^cr^) (q:,/3 = x,y,z), with 
the expectation value taken over an arbitrary quantum 
state of the system. Moreover, invariance oiHxxz under 
the discrete transformations ax — *■ — cr^:, cry — > ctj^, and 
(Jz —Tz implies that G^^„ = and G*/^ = 0. Therefore, 



the element z*^ in Eq 
). The 

spins i and j from Eq. (|3ip yields 



is real, namely. 



1/4 (G^^2; + Gyy). Then, evaluation of the negativity for 



M{L) = — max 



0, dAGl 



^yy 



(33) 

From now on, we will be interested in computing the neg- 
ativity for nearest neighbor spins. The generalized HK 
theorem discussed in Section [nT] implies that we can con- 
sider A as the external potential and (cf crf_|_j^) (for any 
site i) as the relevant density. Thus, J^{L) can be written 
as a function of de/dA for the ground state as well as for 
any minimum energy state in a sector of magnetization 
m (m = 0, ±2, . . . , ±L) and momentum P = (2t:/L)p 
(p = 0, 1, . . . , L — 1). In this direction, it is convenient to 
write the correlation functions G^^^^ in terms of de/dA, 
which results into 



Gl 



i+l _ 



= -2 



dA 



G 



yy 



= -2 £- A 



de 
'dA 



(34) 



A. Ground state entanglement 

For the ground state, we have that Gz = 0. Then, by 
using Eq. ([34]) and Eq. ([33]) . negativity reads 



J\f{L) = -e(i) + (A + 1) 



de{L) _ 1 
dA ~ 2'' 



(35) 



where we have used that |G*'*'''^| < 1 and G. 
G'yl^^ > (Marshall-Peierlsrule). Note that, in 
the energy density e{L) can be seen as a function of 
de/dA by the HK theorem, which is explicitly shown 
in Fig. [H Indeed, this imphes that the negativity can be 
taken as a function of de/dA, which illustrates the dual- 
ity between potential and density established in Eq. ((25l) 
for entanglement measures. The XXZ model is critical 




FIG. 1: (color online) Density energy e as a function of de/dA 
as given by the solution of the Bethe equations for finite size 
chains as well as in the thermodynamic limit (values plotted 
in the range — 1 < A < 1). For finite chains with lattice sizes 
L > 8, the curves get nearly superposed with the curve for 
the infinite chain. Inset: Density energy e versus anisotropy 
parameter A. Note that e can be taken either as a function 
of A (the potential) or de/dA (the density). 

in the interval — 1 < A < 1, with central charge c = 1. 
Then, from Eq. (|35p. we can determine an approximate 
analytical expression for the negativity in terms of en- 
ergy as given by Eq. ([T]). The parameter ^ appearing 
in Eq. ([T|) can be obtained analytically [s^ for the XXZ 
chain, reading 



81117 

7 



(36) 
(37) 



where 7 is defined by 

A = — cos 7 , 7 G [0, tt) 
Then, substitution of Eq. H]) into Eq. (f35|) yields 

l-f A AVTTA 



Mcft{L) = A/'oo + 



TT^C 
67I2 



sm7 - 



7 



(38) 
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where Moo can be computed from Eq. (|35p , with Soo and 
deao/dA directly given by the solution of the model at 
the thermodynamic limit [s^. An exact value for the 
negativity M{L) can be obtained from Eq. ([55)1 by com- 
puting e{L) and de{L) / dA via Bethe ansatz equations for 
each length L. Naturally, this amounts to a much harder 
computational effort for a general A, while Eq. ([55]) di- 
rectly provides the negativity for a finite chain up to or- 
der with no need of solving the Bethe ansatz equa- 
tions for each length L. A comparison between Af{L) 
and Ncft{.L) for 7 = 7r/2 and 7 = 7r/3 is exhibited in 
Tables U and m 



L 


M{L) 


Mcft{L) 


4 


0.457106781187 


0.446378653269 


8 


0.366669830087 


0.366041268056 


16 


0.345995599194 


0.345956921753 


32 


0.340938243195 


0.340935835178 


64 


0.339680713890 


0.339680563534 


128 


0.339366755018 


0.339366745623 


256 


0.339288291732 


0.339288291145 


512 


0.339268677562 


0.339268677525 


1024 


0.339263774123 


0.339263774121 



TABLE I: Comparison between Af{L) and Mcft{L) for 7 = 
7r/2 (the XX model). For an infinite chain, we have that 
negativity is given by 7V(oo) = 0.339262139652. 



L 


N[L) 


Ncft{L) 


4 


0.489830037812 


0.478556230132 


8 


0.401639244141 


0.400889057533 


16 


0.381525197365 


0.381472264383 


32 


0.376621871264 


0.376618066096 


64 


0.375404791436 


0.375404516524 


128 


0.375101148980 


0.375101129131 


256 


0.375025283711 


0.375025282283 


512 


0.375006320673 


0.375006320571 


1024 


0.375001580150 


0.375001580143 



TABLE II: Comparison between Af(L) and Mcft{L) for 7 = 
7r/3. For an infinite chain, we have that negativity is given 
by 7V(oo) = 3/8 = 0.375. 



B. Twisted boundary conditions 

We can also use the results obtained for PBC to in- 
vestigate the finite size corrections to the negativity with 
more general boundary conditions. We will consider here 
the so-called twisted boundary conditions (TBC), which 
can be achieved as the effect of a magnetic flux through 
a spin ring (37| . Remarkably, it has recently been shown 
that TBC may improve multi-party quantum communi- 
cation via spin chains [ssj. In order to consider TBC, 



it is convenient to rewrite the Hamiltonian in Eq. ((29l) 
(with J = 1) in the following form 



H 



xxz 



L 

-2^ 



1 



(39) 

where af = (crj ± icrj')/2 and a^^^ = e±»*crj^ (0 < 
$ < 27r), with $ denoting a phase. The quantum chain 
given by Eq. ([55)1 is solvable by the Bethe ansatz 28] . In 
presence of TBC, Eq. ([T|) still holds, but with an effective 
central c(<I') [285, which is given by 



3$^ 



27r (tt — 7) ' 



(40) 



with 7 defined as in Eq. (1371). Let us take the following 
canonical transformations 13911 



a] (j = l,...,L) 



(41) 



In terms of this new set of operators, the original chain 
with TBC is now given by the periodic chain 



xxz 



1 ^ 

-y 

2 ^ 



(42) 



where cr^_|_i = uf- Defining the operators crj and crj 
through = (ctJ ± icrJ)/2, the Hamiltonian can be put 
in the form 



H 



xxz 



1 ^ 
-Y 

2 ^ 



cos 



(43) 



Note that the Hamiltonian in Eq. (g^l) is both U(l) 
invariant ([^, = 0) and translationally invari- 

ant [Hxxz exhibts PBC in terms of the set {3r^,cr|}). 
Therefore, the two-spin reduced density matrix keeps the 
form given in Eq. pip , with the correlation functions 
replaced by G^^ = ((Ja^^)- Then, the negativity for 
nearest neighbor spins governed by Hamiltonian (I43p can 
be computed similarly as before. By using that Gz = 
(ground state) and |G*'*"'"^| < 1 we obtain 

7V=2max(0, |z|-a), (44) 
where a = {l + Gzz)/4: and z = {G\\+iG±)/A, with Gzz = 

^zz ' '-'ll — '-'xx + ^yy i ^^^^ ^ 1- — xy ~ ^ yx 

(yi). In order to write entanglement in terms of the 
derivatives of the energy density, it is convenient to define 
Hxxz = Hxxz/ coB{<i>/L). Then 



Gil - ?7G_L + AG, 



■)■ 



(45) 
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where e = {H xxz) I L, t] = tan($/L), and A 
A/cos($/i). From Eq. (gl]) we get 

Gp-z — 2 — — , 

dA 



Gi = 2 



de 
dri ' 



^ „ / de -r-de 



(46) 



Therefore, the contribution (|z| — a) for expression for the 
negativity in Eq. reads 



\z\ — a = — 




de de 
dA 



(47) 



In order to obtain the results in terms of $ and A, we 
make use of the expressions 



de 
drj 
de_ 



eos(| 

de_ 
dA' 



, de /$ 



(48) 
(49) 



Hence, finite size corrections to entanglement can be 
found now by using Eq. ([T]) [replacing the central charge 
c by the effective central charge c($) as in Eq. (j40|) ] into 
Eq. (|T7)) . Examples comparing the negativity Mcft{L) 
for nearest neighbors up to o{L~^) and the exact value of 
the negativity JV{L) (obtained through the numerical so- 
lution of the Bethe ansatz equations derived in Ref. [28l |) 
are exhibited in Tables [1111 and HVl below. 



L 


AA(L) 


AfcFT{L) 


4 


0.406774810601 


0.446378653269 


8 


0.354315234931 


0.366041268056 


16 


0.342922395530 


0.345956921753 


32 


0.340170924101 


0.340935835178 


64 


0.339488945731 


0.339680563534 


128 


0.339318816833 


0.339366745623 


256 


0.339276307427 


0.339288291145 


512 


0.339265681501 


0.339268677525 


1024 


0.339263025109 


0.339263774121 



TABLE III: Comparison between the exact evaluation of 
N[L) and the approximate expression Mcft{Ij) (up to or- 
der L"^)) for 7 = 7r/2 and 3> = 7r/2. 

Note from TablesHlandHTTlthat. for 7 7r/2, Ncft{L) 
gives the same result either for $ = (PBC) or $ = 7r/2, 
which is an indication that TEC should not affect the 
negativity (up to o{L~^)) in the case of the XX model. 
Indeed, this can be analytically proved. In this case, the 



L 


U{L) 


Ucft{L) 


4 


0.400000000000 


0.452230707893 


8 


0.381121448251 


0.394307676973 


16 


0.376577662094 


0.379826919243 


32 


0.375405200439 


0.376206729811 


64 


0.375102994373 


0.375301682453 


128 


0.375025983614 


0.375075420613 


256 


0.375006526789 


0.375018855153 


512 


0.375001635654 


0.375004713788 


1024 


0.375000409414 


0.375001178447 



TABLE IV: Comparison between the exact evaluation of 
Af{L) and the approximate expression Mcft{L) (up to or- 
der L"^)) for for 7 = 7r/3 and $ = 27r/3. 



anisotropy is A = 0, which implies that 7 = 7r/2 and 
^ = 2. Then, from Eqs. (P) and (gOl), we have 



e{L) 


— ^00 


de 


2 $ 


d^ 




de 


deoo 


dA 


dA 



A=0 



2 
3i2 



+ — + o(L-2). 



(50) 



By inserting the above equations into Eqs. (gHl) and (^5]) . 
it can be shown that the negativity as given by Eq. ((44)) 
gets 



cft 



where |e 



|£($=o) I + 



de„ 



OA 



A=0 



■i,-l + oiL-% (51) 



(i>=o) 



1 = 



2eoo7rC/(6L2). Hence, Eq. §1 



implies that the negativity for the XX model with TBC 
is not affected by the phase $ up to order l/L^. 



C. Excited states 

Let us consider now the structure of the negativity 
for the excited states in the XXZ model with PBC. The 
U{1) and translation symmetries allow us the decompo- 
sition of the associated eigenspace of Hxxz into dis- 
joint sectors (fixed magnetization and momentum) la- 
belled by the quantum numbers r = 0, 1,2, and 
p=l,2,...,L — 1, which give the number of spins up 
in the cr^ basis and the eigenvalue of the momentum 
P ~ {2n/L)p, respectively. An exact evaluation of the 
negativity for nearest neighbor spins can be performed 
from Eq. (|33p by taking a non-vanishing magnetization 
density Gz and by using Eq. (|34p . where the energy of 
the excited state is obtained through the solution of the 
Bethe ansatz equations. This is illustrated in Fig. [21 
where we plot the negativity between nearest neighbors 
in a chain of length L = 256 sites for minimum energy 
states with zero momentum in several magnetization sec- 
tors. These states have anomalous dimensions Xn given 
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0.3 



N0.2 



0.1 



n = 
n = 4 
n = 8 



I ' I ' I ' I ' I 








-0.4 -0.2 0.2 0.4 



-0.5 





A 



0.5 



FIG. 2: (color online) Negativity for minimum energy states 
as a function of the anisotropy A for L = 256 sites. Note 
that the curves are nearly the same, indicating a unique en- 
tanglement pattern in the critical region. Inset: Negativity 
as a function of A in a larger zoom scale. 



by m 



: (tt - 7) 
2^ ' 



(52) 



where n = L/2 — r and j + j' = p = Q in Eq. Re- 
markably, note that the negativities for the minimum 
energy states plotted are nearly the same, indicating a 
unique entanglement pattern in the critical region. In- 
deed, this is a more general result, which holds also for 
other excited states. For instance, let us take the so- 
called marginal state [l^, which is a state that will be 
taken in the sector n = with anomalous dimension 
X = 2 (independently of 7) and = 0. Exact compu- 
tation in Table |V] below shows that its negativity is also 
close to the values found in Fig. O Indeed, we can show 



A 


Marginal State 


Ground State 


0.505023772863 


0.265447369819 


0.266151418398 


0.205023772863 


0.315358910123 


0.316005319520 


0.005023772863 


0.338660739066 


0.339288291732 


-0.204976227137 


0.357090720201 


0.357706303640 


-0.504976227137 


0.374867541783 


0.375473489099 



TABLE V: Comparison between A/'(A) for the marginal state 
and and A/'(A) for the ground state in a chain with L = 256 
sites. 



at order L"^ [see Eqs. ^ and ^] . According to Eq. p3|) . 
such a difference is also reflected in the negativity of near- 
est neighbor spins, which explains the behavior displayed 
in both Fig. [2] and Table fVl This can explicitly be shown 
by inserting Eq. ([2]) into Eq. ((33|) . As an illustration, we 
take the minimum energy states with zero momentum in 
a given magnetization sector labelled by n. For this case, 
the negativity can be evaluated as 



AA„(L)=AC + 



IG. 



67L2 



3771^ -I- TT sin 7 Goo 



-7r(AGoo + |Goo|) ( — - 2nC0t7 



o{L- 



where 



7r(c-6n2), = |Goc| 

(tTC - 6n^(7r - 7)) , Goo = 



■de'^/dA-1/2, 
+ Ade^/dA. 



with ejf denoting the energy density of the excited state 
as L ^ 00. Note that this unique pattern of entangle- 
ment, which has been explicitly derived here, is in agree- 
ment with the general discussion of Sec. IIVI This is in- 
deed exhibited in Eq. (I28p . Naturally, similar expressions 
can be obtained for excited states higher than the mini- 
mum energy states. 



VI. CONCLUSION 

In conclusion, we have investigated the computation of 
finite size corrections to entanglement in quantum criti- 
cal systems. These corrections were shown to depend on 
the central charge of the model as well as the anomalous 
dimensions of the primary operators of the theory. Our 
approach has naturally arisen as a general consequence 
of the application of CFT and DFT methods in critical 
theories. This framework has been illustrated in the XXZ 
model, where we have shown that: (i) entanglement in 
spin chains with arbitrary finite sizes can be analytically 
computed up to order o{L~'^) with no need of solving the 
Bethe ansatz equations for each length L; (ii) Conformal 
towers of excited states displays a unique pattern of en- 
tanglement in the critical region. Indeed, we have been 
able to provide a general argument according to which 
this unique pattern of entanglement should appear in 
all conformally invariant models. Further examples in 
higher dimensional lattices and higher spin systems are 
left for future investigation. 



that entanglement in the critical region of the XXZ chain 
will exhibit a unique pattern for all states accessible via 
the CFT associated with the model. As discussed in Sec- 
tion ini each primary operator of the theory corresponds 
to a tower of states with energies given by Eq. All 
these states in the towers will have energies which differ 
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